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Abstract 

A class of non-linear eigenvalue problems defined in the form of operator poly- 
nomials is investigated. The problems are related to wave equations which appear 
in a relativistic quantum field theory. Spectral asymptotics for this class are found 
explicitly. The properties of operator polynomials are analyzed for scalar, spinor 
and gauge fields. It is also shown how to use these results in finite temperature 
theories. 
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1 Introduction 



Consider a class of eigenvalue problems which have the following form: 

P(u)<j> u EE (P n + UJP^ + ... + U n ~ X P X + 0J n P )^ = , (1.1) 

where ui is a complex spectral parameter and 0^ is an eigenfunction. Pk (k — 0, n) are 
partial differential operators of fc-th order acting on a Hilbert space H. P{oo) represents 
an example of what is called the operator polynomial or polynomial operator pencil [1]. 
In general Pk and Pi do not commute if k ^ I and (1.1) is essentially non-linear problem 
with respect to the eigenvalue u>. 

Problems like (1.1) appear in different applications (oscillations in a viscous fluid, 
Schrodinger equation with energy dependent potential and etc), see [1] for the introduction 
and references. The focus of the present work is on applications related to quantum field 
theory. 

Consider equation of motion for a free field <p(x, x°) (x and x° denote spatial and time 
coordinates, respectively). It can be always written as 

P(id )<f>(x,x°) ee (P 2 + P^ido) + P o (^o) 2 )0(x,x°) = , (1.2) 

where d = d/d . Operators Pk depend only on spatial derivatives and may have ad- 
ditional indexes related to a spin-tensor structure. Field <fi propagates in an external 
classical background. If the background is stationary the variables can be separated by 
choosing <f>(x,x°) = e~ tuJX ° <f> u (x) . From (1.2) one gets 

P{u)<f> u {x) = (P 2 + uP-y + ^Po^ix) = . (1.3) 

This is a quadratic operator polynomial, a particular case of problem (1.1). 

Eigenvalues ui determine the spectrum of energies of single-particle excitations which 
is used to compute a number of important physical quantities, like the vacuum (Casimir) 
energy E — ±| M, free energy and etc. 

Consider for instance a charged Dirac field in an external static electromagnetic 
field with the potential = A^x). The equation is 

[7^^ + 771^ = , (1.4) 

where = V M — ieA^ and e is the charge. To get for the single-particle spectrum 
problem like (1.3) one has to act by the operator [7 M -D M — m] on the left hand side (l.h.s.) 
of (1.4) and make substitution ip(x,x°) = e~ twx ip u (x). This yields 

[D k D k -m 2 + e 2 A 2 - e ~[T - 2t 1 ° 1 k d k A ) + uj 2eA + uj 2 ]^{x) = , (1.5) 

where T = (djAk — dkAj)^^ k and k,j = 1, 2, 3. When the time component Aq of the po- 
tential is a function of spatial coordinates (1.5) cannot be reduced to a standard eigenvalue 
problem because operator P 1 = 2eA does not commute with P 2 . 
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Sec. 2 of this contribution contains a summary of results obtained in [2, 3]. The aim 
is to show that the notion of spectral geometry can be extended to a class of quadratic 
operator polynomials which appear in problems of quantum field theory. Applications to 
finite temperature theories are discussed in Sec. 3 where, as an example, the obtained 
results are used to describe the Debye screening in hot plasma. Finally, some important 
features of operator polynomials for scalar, spinor and gauge fields are analyzed in Sec. 
4. 



2 Spectral Geometry of Quadratic Polynomials 

To describe the results we rewrite (1.3) in the form 

W 2 -L( W )]i(x) = , (2.1) 

L(u) = Lz + uLi+^Lo , (2.2) 

where L(u>) = uj 2 — P{oj) and L k are partial differential operators of the fc-th order. 
Suppose for a moment that w is a complex parameter and consider a one-parameter 
family of operators L{uj). Let A k {uj) be eigenvalues of L(u>) for a given uj, 

LM^=A*( W )tfff. (2.3) 

If the spectrum A k (oj) is known for any uj, the spectrum of problem (2.1) is determined 
by roots of equation 

X (u,A k )=0, (2.4) 
x(w,Ajfc) = uj 2 - A k (u). (2.5) 
Denote these roots by Uk,% and define the following function on the spectrum 

x'Ki) = ^xKi,Aib) where u 2 ki = A fc (w fc)i ). (2.6) 

It is assumed that at first the derivative in x'i^k^) over uj is taken for the fixed branch of 
eigenvalues A k {uj) and after that the result is considered on the corresponding root u k ^. 
In what follows we write uo instead of uj k ^. 

We make three assumptions: i) L{uj) is a Laplace type operator of the form: 

L(u) = -(V fc + iA k + iua k )(V k + iA k + iua k ) +uB + V (2.7) 

which acts on sections to vector bundles over a <i-dimensional compact Riemannian man- 
ifold Aid', ii) the spectrum of L(0) is strictly positive; iii) the function x\ tjJ ), see (2-6), is 
positive (negative) on positive (negative) eigenvalues uj. 

First two conditions are technical (see details in [3]). The origin of condition (iii) is 
explained in Sec. 4. At this point we note that in the limit when the polynomial is trivial, 
L{uj) = L 2 , this condition is satisfied because x'i 00 ) = 2^- 
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To describe the spectrum of (2.1) at large we introduce the pseudo-trace 

K(t)^We-^\ t>0, (2.8) 

where summation goes over the real eigenvalues cu. The asymptotics of K(t) at small t is 
related to the distribution of large Given conditions (i)-(iii) the following asymptotic 
series takes place at small t [3] 

1 oo 

K V~J^f o [ a " tn + b « tn+1/2 }- ( 2 - 9 ) 
a n and b n can be computed by using expansion formula for the heat kernel 

1 oo 

K u (t) = Tr e -«*"> ~ £ [o^t* + W* n+1/2 ] (2.10) 

which is well-known [4]. For operators defined in (2.7), 

n n 
a n(uj) = a m,nU m , b n (u) = ^ b m ^ m (2.11) 

m=0 m=0 

and it can be shown that [3] 

2» n _ m r(-f + m) 

a « = XI d r a 2(m-n),m, (2.12) 

m=n 1 ( — 2 H ) 

b n= m / d _! rh(m-n),m- (2.13) 

m=n 1 ( 2 \~ Tl) 

Coefficient a is proportional to the volume of M d , other a n , 6 n are integrals of local 
functionals of background fields and b n are non-trivial when M d has a boundary. The 
notion of spectral geometry does extend to a class of quadratic operator polynomials. 
The proof of (2.9), (2.12), (2.13) is based on relation [3] 

K{t) = I" 'duo f dz e -<" 2 (*-«) (2u + -dj) (K u (iz) + K^iz)), (2.14) 

where a contour C goes in the complex plane just below the real axis Im z = 0. Thus K(t) 
is determined by part of K w [t) which is symmetric with respect to the transformation u 
to -co. This property is a direct consequence of condition (iii) which we imposed on the 
spectrum of the operator polynomial. This also explains why coefficients a n and b n in 
(2.12), (2.13) are determined by the symmetric part of coefficients a n (u) and b n {uj). 
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3 Applications to Finite- Temperature Theories 



Spectral asymptotics of quadratic operator polynomials can be used in a number of phys- 
ical applications. One of them is studying the behaviour of the free energy of non- 
interacting quanta at high temperatures. The free energy for Fermi particles at temper- 
ature T = (3~ l is defined as 

ns) = -/r 1 j> (1 + e-^*) . t 3 - 1 ) 

k 

Summation goes over single-particle energies E k = |u; fe | which are determined by the real 
eigenvalues uj k of the corresponding wave equations. 

It is convenient to start with the case when the single-particle spectrum is determined 
by the standard linear eigenvalue problem, i.e., when operator L{uj) = L 2 and it does not 
depend on uj. This case is well studied. It is known that in three dimensions d = 3 the 
free energy in the high temperature limit has the following asymptotic expansion [5] 

m A) ~ ~ 7 ^a - ^TV - JL ln(7»a 2 + 0(T~ 2 ) (3.2) 

(we neglect boundary terms for simplicity). Here a n are the coefficients of the asymptotic 
expansion of the trace Tr e~ tL2 , g is a dimensional parameter. 

In the high temperature limit F(/3) is determined by behaviour of the single-particle 
spectrum at large \u\. For polynomial (2.1) this behaviour is described by asymptotic 
formula (2.9) which has exactly the same form as expansion of Tr e~ tL2 . Therefore, 
for systems whose spectrum is determined by (2.1) formula (3.2) can be used if a n are 
calculated with the help of (2.12) (analogously for boundary terms if they are present). 

Let us use this method to derive the free energy of electron-positron plasma in an 
external static electromagnetic field A^x). (Finite-temperature theories in gravitational 
backgrounds are discussed in [2]). Coefficients a n are determined by (2.12) for problem 
(1.5). The corresponding operator is 

L(u) = -V k V k + m 2 - e 2 A 2 + - 2i 1 ° 1 k d k A ) - 2eA Q u . (3.3) 
By using standard results [4] one easily finds coefficients of L(u) and obtains 

= ao,o = 4 / d 3 x , (3.4) 

Jm 3 

ai = a ,i + \a 2 ,2 = 4 / d 3 x{-m 2 + 2e 2 A 2 ) . (3.5) 

2 J M 3 

1 3 2e 2 r 

a 2 = a ,2 - -a 2 , 3 + 704,4 = — / d^xF^F^ , (3.6) 

2 4 3 Jm 3 

F^ u = d^A u — d v A^. Consider the effective thermodynamic potential of the system at 
high temperature by neglecting the zero-point energy 

Q((3, A) = — L(A) + F((3, A) = \ d 3 xF^ + F((3, A) , (3.7) 
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where L(A) is the classical Lagrangian of the background potential A^. Eqs. (3.2), (3.4)- 
(3.6) yield 

m A) = - j M3 d 3 x {-^-F^ + \m\T)A^ , (3.8) 

c(T) = l-^ln(T/p) , m 2 (T) = ie 2 T 2 . (3.9) 

These formulas reproduce the result which is well-known in plasma physics. The last 
term in (3.8) is responsible for screening effects in hot plasma and m(T) is the correct 
expression for the Debye mass. 

Note that conventional derivation of (3.8) is based on computation of one-loop dia- 
grams at finite temperature in the limit of vanishing gauge fields. Our calculation is done 
for arbitrary static gauge potential. 

Coefficient a 2 in (3.6) has a Lorentz invariant structure. For this reason the logarithmic 
term in expansion (3.2) results in a finite renormalization of the background potential 
(see the first term in r.h.s. of (3.8)). The Lorentz invariant structure of a 2 is not a mere 
coincidence. It reflects a more general property established in [3]: coefficient a 2 of the 
pseudo-trace expansion coincides with the coefficient A 2 of the asymptotic expansion of 
the heat kernel of the covariant four-dimensional operator P(id ) = — <9q — L(id ). 

4 Operator Polynomials for Different Field Theories 

Spectral asymptotics given in Sec. 2 are obtained by using certain conditions. In partic- 
ular, the structure of coefficients a n and b n in (2.12), (2.13) is determined by condition 
(iii). Let us analyze the origin of this condition for models of non- interacting fields. 

Consider a pair / 1; f 2 of solutions to field Eq. (1.2) (or (1.4), for Dirac case) and 
introduce the relativistic product, a Hermitian bilinear form (fi, f 2 ) = (fi,fi)* defined 
on a space-like hypersurface E in a hyperbolic spacetime M. as 

(/i,/ 2 > = Jj^j,(f 1 ,f 2 ) • (4.1) 

The "current" j M (/i, f 2 ) is divergence free, = 0, and (4.1) does not depend on E. 

If the background is stationary one can introduce a conserved canonical energy H[f] 
for a given solution /. It can be shown that for free scalar, spinor and gauge fields the 
energy can be written in the form 

H[f]= l -(f,d f)+c.c, (4.2) 

where "c.c." is a term obtained from the first term in the r.h.s. of (4.2) by complex 
conjugation. It follows from (4.2) that for a solution f u with certain frequency cu (idof^ = 

H\f u ]=u(f u ,f u ) . (4.3) 
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In general the system may have excitations with complex frequencies uj. It is easy to 
see that they have vanishing relativistic norm (f u , fj) and, hence, the vanishing energy. 
Such modes along with zero-frequency modes are not quantized and should be considered 
separately from real-frequency solutions. This is the reason why pseudo-trace (2.8) is 
defined for real uj. 

4.1 Scalar Fields 

Consider integer spin fields. To have a well-defined theory we require that canonical 
energy is non- negative, 

H[f u ]>0 , (4.4) 

for any uj. If there are single-particle excitations with negative energy H the Fock space 
in the corresponding quantized theory contains states with negative energies unbounded 
from below. 

Condition (iii) of Sec. 2 for scalar fields <fi follows from (4.4). Indeed it is possible to 
show [3] that in this case 

=X / H(^ w ,0 w ) , (4-5) 

where x'O^) was introduced in (2.6) and (<f> u ,<f> u ) = J dVlcp^l 2 is a positive-definite norm 
in a Hilbert space of single-particle wave functions. (4.4) together with relations (4.3), 
(4.5) implies that x'i 00 ) = £ ( UJ )\x'{ UJ )\ where e{uj) is the sign function. 

4.2 Gauge Fields 

We discuss now theory of non-Abelian gauge fields with the gauge group SU (N). Let _B M 
be a solution to Yang-Mills equations, D^Fffl = 0, where Fffl is the strength tensor of 
(we work in the adjoint representation). In the linear order, perturbations of the 
Yang-Mills field on the background obey the equation 

V U V U A^ - V^Av + 2iF^A v = , (4.6) 

where X> M = + iB^. It is convenient to work in the background gauge V^A^ = where 
equation transforms to 

Trv„A lt + 2iFWA v = . (4.7) 

Fixing this gauge leaves a freedom in the gauge transformations 5A^ = where A is a 
solution to 

V li W i \ = Q . (4.8) 

Suppose now that B^ is a stationary field. Then Eqs. (4.6)-(4.8) result in operator 
polynomials of the form (1.3). Let u be the physical spectrum of (4.6), the spectrum 
related to vector Eq. (4.7) and U( ) the spectrum of (4.8). It is clear that u;'s represent 
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a subset among eigenvalues U(iy The pseudo-trace for the physical spectrum can be 
represented as 

K(t) = ^e-^=K w (t)-2K i0) (t)= £ e~ t "h-2 £ e"S» ■ (4.9) 

(For SU(N) group Eqs. (4.6)-(4.8) are invariant with respect to the charge conjugation 
and their spectrum is symmetric with respect to change u to — uj.) The double subtraction 
of the sum over u;(o) in (4.9) has the following explanation. One subtraction eliminates 
solutions of (4.7) which do not respect the gauge conditions, V^N 1 = 0^0. It is easy to 
see that such modes have frequencies because obeys scalar Eq. (4.8). Additional 
subtraction of K^(t) in (4.9) eliminates pure gauge solutions = V^X. The last term 
in r.h.s of (4.9) can be interpreted as contribution of ghosts. 

There is an analog of formula (4.5) for transverse solutions (V^A^ = 0) to (4.6). The 
relativistic norm for the solutions with energy uj can be written as 

(AmA^HxVXAmA,) , (4.io) 

where (A W1 A U ) = J dV(A+) fl (A u ) fl and quantity x'i 00 ) is defined for the polynomial cor- 
responding to vector problem (4.7). The residual gauge freedom can be used to impose 
condition Aq = 0. For such solutions the norm [A^^AJ) becomes positive-definite. 

Suppose that the canonical energy of perturbations A^ on the background is non- 
negative, iJyM^4] > 0. Require also the same property for canonical energy of the 
scalar field described by (4.8), H[X] > 0. Then one can show that operator polynomials 
associated with vector, (4.7), and scalar, (4.8), problems obey condition (iii) of Sec. 2. 
The proof is based on relations (4.5), (4.10). 

4.3 Dirac Field 

As is known, the classical energy for the Dirac field is not positive-definite. Thus, the proof 
of condition (iii) in this case should be different. Consider the Dirac field in Minkowsky 
space-time in the presence of a static gauge potential A^. The relativistic norm of a 
single-particle wave function is = J dV if) + ip. If the field has a mass, m/ 0, one 

can prove the relation [3] 

where (V'w,^) — / dV ip^i^oip^ and 70 is ant i- Hermit ian. It follows from (4.11) and 
positivity of the norm that ^{uj) does not vanish for massive fields. Now, if 

A^ = 0, one has x'i 00 ) = 2oo and condition (iii) is satisfied. The validity of (iii) for A^ ^ 
can be proved by continuity. Because x'i 00 ) quantity x'i 00 ) cannot change the sign 
when A a is switched on from zero to some value. 
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